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A  number  of  mechanisms  are  known  to  be  responsible  for  the  inter¬ 
dependence  of  stratification  and  shear.  These  include  interleaving, 
internal  waves,  and  instabilities  such  as  double  diffusion  and 
breaking  internal  waves.  In  the  strongly  stratified  seasonal 
thermocline  the  shear  behaves  approximately  as  S2®N2  (Grabowski,  1980) . 
Patterson  et.  al.  (1981)  showed  that  N2  and  S2  are  better  correlated 
over  vertical  scales  of  30  meters  or  more.  This  calculation  of  cross 
coherence  between  simultaneous  profiles  of  N2  and  S2  was  intended  to 
investigate  the  possibility  of  a  limiting  vertical  length  scale  below 
which  the  correlation  of  N2  and  S2  is  small  or  zero,  and  above  which 
the  correlation  is  good.  On  the  basis  of  a  limited  data  set  (9  YVETTE 
rofiles) ,  N2  and  S2  appear  in  general  to  be  well  correlated  at  wave- 
engths  larger  than  about  5-10  meters,  except  in  the  presence  of  large 
vertical  geostrophic  shear.  These  results  offer  some  further  evidence 
that  levels  of  shear  activity  may  be  inferred  from  stratification  at 
scales  down  to  about  10  meters 
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Section  1 
INTRODUCTION 

It  has  been  well  documented  in  recent  years  that 
stratification  and  shear  in  the  ocean  are  often  highly  inter¬ 
dependent.  A  number  of  mechanisms  are  known  to  be  respon¬ 
sible  for  this  interdependence,  including  interleaving, 
internal  waves,  and  instabilities  such  as  double  diffusion 
and  breaking  internal  waves.  For  example,  Eriksen  (1978) 
documented  the  effects  of  breaking  internal  waves  on  the 
relationship  between  shear  squared  (S2)  and  Vaisala  frequency 
squared  (N2).  His  scatter  plots  of  N2  and  S2  computed  over  a 
vertical  interval  of  6.3  m  show  the  limiting  Richardson 
number  Ri  *  n2/S2  to  be  about  1/4.  Shear  instability,  there¬ 
fore,  seems  to  bound  S2  to  values  less  than  the  local  value 
of  4N2. 

Eriksen' s  observations  show  that  Richardson  number 
approaches  1/4  rather  infrequently.  Internal  waves  are  a 
possible  mechanism  for  the  apparent  correlation  between  S2 
and  N2  when  Ri  >  1/4.  Johnson  and  Sanford  (1980)  showed  that 
an  anisotropic  internal  wave  field  superimposed  on  a  vertical 
temperature  gradient  resulted  in  significant  coherence 
between  vertical  shear  and  temperature  gradient.  WKB  scaling 
arguments  suggest  that,  at  least  in  the  deep  ocean,  the  shear 
due  to  linear  internal  waves  behaves  as 

S2  «  N3. 

Observations  made  by  the  free-fall  profiler  YVETTE  below  the 
seasonal  thermocline  tend  to  support  this  proportionality. 
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In  the  strongly  stratified  seasonal  thermocline,  the  shear 
behaves  more  nearly  as 

S2  «  N2 
(Grabowski,  1980). 

The  similarity  between  profiles  of  N2  and  S2  was 
first  pointed  out  by  Simpson  (1975).  Figure  1.1  shows 
profiles  of  N2,  S2,  and  Ri  from  the  free-fall  shear  profiler 
YVETTE.  In  the  seasonal  thermocline  (between  the  dashed 

lines),  and  to  a  lesser  extent  below  the  seasonal  thermo¬ 
cline,  Ri  is  limited  by  the  value  1/4.  The  similarity 
between  the  profiles  of  N 2  and  S2  is  immediately  obvious. 

In  Patterson  et  al_.  (1981),  we  presented  linear 
cross-correlation  coefficients  computed  between  N2  and  S2, 
for  different  portions  of  an  YVETTE  profile.  We  showed  that 
smoothing  the  profiles  of  N2  and  S2  generally  improves  corre¬ 
lations.  The  implication  is  that  N2  and  S2  are  better  corre¬ 
lated  over  long  vertical  length  scales  than  over  short 

scales.  This  calculation  led  to  the  question  of  whether 

there  is  a  limiting  length  scale,  above  which  N2  and  S2  are 

well  correlated,  and  below  which  the  correlation  is  poor.  We 
anticipated  that  a  calculation  of  coherence  between  N2  and  S2 
would  provide  a  mo’re  definitive  answer  than  simply  varying 
the  smoothing  of  the  cross  correlation.  Thus  we  tried  the 

coherence  calculation  on  nine  different  YVETTE  profiles. 
This  technical  note  describes  the  results  of  that  brief 
investigation. 

In  Section  2  we  present  an  outline  of  the  method 
used  to  compute  coherence.  In  Section  3  we  present  the 

results  and  a  brief  discussion. 
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N2  -  S2  COHERENCE:  METHOD  OF  COMPUTATION 
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The  locations  of  the  nine  stations  for  which  coher¬ 
ence  spectra  were  computed  are  listed  in  Table  2.1  and  shown 
in  Figure  2.1.  The  data  used  were  from  the  depth  intervals 
below  the  mixing  layer,  and  included  both  the  strongly  strat¬ 
ified  seasonal  thermocline  and  the  more  weakly  stratified 
layers  below.  The  data  were  sampled  at  approximately  1  meter 
intervals.  Sixteen  meters  was  chosen  as  the  longest  wave¬ 
length  of  interest,  and  thus  each  station  was  divided  into 
records  of  16  samples  (approximately  16  meters)  long.  The 
number  of  non-overlappinq  16  sample  records  in  a  given 
station  is  denoted  by  the  integer  L. 

We  next  outline  the  method  used  for  computing 

coherence  between  N2  and  S2.  Data  points  in  the  time  series 

2  2 

of  N2  and  S2  are  denoted  Nj  and  Sj  respectively. 

1.  Truncate  the  original  station  data  records  to  a 
length  M,  such  that 

M  =  LJ, 

where  J  is  a  power  of  two,  and  L  is  the  integer 
defined  above.  We  then  have  data  records  Nj 
and  Sj  where  j  =  1,2,...,M. 

2.  Subtract  the  mean  from  each  of  these  records. 

3.  Divide  the  station  data  record  into  2L  -  1 

overlapping  intervals.  Apply  a  cosine  taper 
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Table  2.1 
YVETTE  STATIONS1 


Station 

Number 

Time 

(GMT) 

Date 

Latitude 

(N) 

Longitude 

(W) 

Goranent 

5 

1942 

5  Nov.  75 

32619' 

64*34' 

Near  Bermuda 

8 

0225 

8  Nov.  75 

35°00' 

66*30' 

Sargasso  Sea 

9 

1219 

li 

11 

11 

it 

10 

1814 

9  Nov.  75 

38  *09 ' 

69*06' 

Gulf  Stream 

11 

0036 

10  Nov.  75 

38°05' 

69*03' 

12 

1312 

11 

38*15' 

69*07’ 

11 

18 

— 

7  May  77 

22*47' 

70*43' 

Edge  of  thenrocline  eddy 

21 

— 

9  May  77 

22°27' 

70*57' 

Center  of  thermocline  eddy 

23 

— 

16  May  77 

36°24'  ■ 

67*36' 

Outer  part  of  GSR2 

» 


1  Adapted  from  Lambert  _et  al.  (1980) 

2  Gulf  Stream  Ring 


» 


I 


2-2 


Figure  2.1.  Positions  of  YVETTE  stations 
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(Hanning)  window  to  the  data  in  each  of  these 
intervals. 

4.  Compute  the  complex  valued  Fast  Fourier  Trans¬ 
forms  : 


for  each-  of  the  overlapping  intervals.  Here 

N~  and  denote  spectral  coefficients, 
n  n 


5.  Compute  the  autospectra  (G^,  Ggi ,  cospectra 
(P)  and  quadrature  spectra  (Q): 


2h 

J 

tN 

Gs(£n)  ' 

2h 

J 

i^r 

P(t„)  = 

2h 

(Re^n 

ReS„  + 

imN^lmS^  ) 

J 

n 

n  n  f 

Q(fn)  - 

2h 

J 

(R^ 

“*£  - 

ImN^ReS~)  , 
n  n ' 

where  h  is  the  sampling  interval  and  fR  « 
for  n  —  0,1,2,... , -j*  1  . 
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6.  Compute  the  coherence  amplitude  and  phase: 


f 


"<p(fn)>2  +  <Q(fn)>2"j 

<W>  J 


where  the  angled  brackets  indicate  averages 
over  the  2L-1  autospectra,  cospectra,  or  quad¬ 
rature  spectra. 


When  coherence  is  computed,  an  implicit 
assumption  is  made,  namely,  that  the  relation¬ 
ship  N2  and  S2  is  linear.  The  coherence 
amplitude  is  a  measure  of  the  correlation 
between  N2  and  s2  in  each  frequency  band. 
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Section  3 

RESULTS  AND  DISCUSSION 


^  t 


The  coherence  and  phase  spectra  for  the  nine 
stations  are  shown  in  Figures  3.1  -  3.9.  Each  coherence 

spectrum  was  computed  with  L  ~  30  indicating  approximately  60 
degrees  of  freedom.  The  level  of  significance  at  90  percent 
confidence,  determined  according  to  Koopmans  (1974),  is 
denoted  by  a  dashed  line  in  each  figure.  The  point  at  0  cpm 
can  be  interpreted  as  the  coherence  at  wavelengths  of  16  m 
(the  total  record  length).  The  range  covered  by  this 
analysis  is  thus  about  16  m  to  2  m. 

For  all  but  Station  12,  the  phase  is  stable  and 
near  zero  in  the  region  of  high  coherence  confirming  that 
spatial  fluctuations  in  N2  and  S2  are  synchronized  at  scales 
longer  than  4-10  meters.  This  is  consistent  with  the 
improvement  in  cross-correlation  between  N2  and  S2  profiles 
reported  by  Fatterson  et.  al.  (1980)  when  the  profiles  were 
first  low  pass  filtered  with  a  triangular  weighting  function 
of  5  meter  half-width  (approximate  wavelength  cutoff  of  10 
meters) . 

The  results  are  not  unambiguous.  Station  5  and  8 
show  very  high  coherence  (~  0.7  -  0.8)  at  small  wavenumbers 
with  a  sharp  drop  at  wavenumbers  above  about  0.2  cpm. 
Station  9  shows  significant  coherence  (but  with  an  amplitude 
of  only  about  0.4)  and  with  a  drop  at  wavenumbers  above  0.1 
cpm.  Stations  5,  8,  and  9  were  all  made  in  almost  the  same 
location  in  the  central  Sargasso  Sea.  Station  10  (in  the 
Gulf  Stream)  closely  resembles  Station  9.  However,  Stations 
11  and  12,  also  from  the  Gulf  Stream,  show  insignificant 

3-1 
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Coherence  naplttude  and  phas t  spectra  for  YVETTE  Station  S.  Thr 
dashed  line  Is  the  level  of  significance  at  905  confidence  for 
coherence  amplitude.  955  confidence  Hefts  on  phase  are  •  9°. 


but  for  YVETTE  station  10 


station  IS 


coherence  at  small  wavenumbers  with  significant  coherence 
only  between  about  0.1  and  0.2  cpm,  in  sharp  contrast  with 
Station  10. 

Stations  18  and  21  are  from  the  same  location,  near 
the  Bahamas.  Both  exhibit  similar  coherence  properties, 
namely,  moderately  high  coherence  (0.65)  at  small  wave- 
numbers,  with  sharp  drops  in  coherence  below  0.1  cpm  (Station 
18)  and  0.18  cpm  (Station  21).  The  coherence  values  from 
Station  23  (near  the  edge  of  a  Gulf  Stream  Ring)  at  small 
wavenumbers  are  not  significant. 

There  is  an  overall  trend  that  suggests  that  small 
wavenumber  coherence  is  higher  in  regions  of  relatively  low 
geostrophic  shear.  However,  the  wavenumber  above  which  the 
coherence  becomes  insignificant  is  not  .well-determined  in  any 
case.  The  best  that  can  be  said  is  that  it  occurs  somewhere 
between  0.1  and  0.25  cpm. 

In  summary,  on  the  basis  of  analysis  of  a  very  few 
profiles,  N2  and  S2  appear  in  general  to  be  well  correlated 
at  wavenumbers  smaller  than  about  0. 1-0.2  cpm,  except  in  the 
presence  of  large  vertical  geostrophic  shear.  Profiles 
obtained  in  regions  of  large  vertical  geostrophic  shear  show 
no  significant  coherence  at  the  smallest  wavenumbers 
addressed  in  this  analysis  (~  0.06  cpm).  In  contrast,  the 
results  of  Patterson  et.  al_,  (  1980)  suggest  that  such 
profiles  (e.g.,  Station  12)  are  well  correlated  at  longer 
wavelengths  (smaller  wavenumbers).  One  station  (Station  10), 
obtained  in  the  Gulf  Stream,  does  exhibit  coherence,  in  the 
present  analysis,  at  wavelengths  of  16  meters  and  is  thus  an 
exception.  These  results  do  offer  some  further  evidence  that 
levels  of  shear  activity  may  be  inferred  from  stratification 
at  scales  down  to  about  10  meters. 
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